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Abstract 

Generalizations of the Kovalevskaya, Chaplygin, Goryachev-Chaplygin 
and Bogoyavlensky systems on a bundle are considered in this paper. More- 
over, a method of introduction of separating variables and action-angle vari- 
ables is described. Another integration method for the Kovalevskaya top on 
the bundle is found. This method uses a coordinate transformation that re- 
duces the Kovalevskaya system to the Neumann system. The Kolosov analogy 
is considered. A generalization of a recent Gaffet system to the bundle of Pois- 
son brackets is obtained at the end of the paper. 

Dedicated to the 150-th anniversary of S. V. Kovalevskaya 

It was the 150-th anniversary on January 15, 2000 of Sofia Vasilievna Ko- 
valevskaya, an outstanding women-mathematician, whose works made a great con- 
tribution into development of the world science. 

Her literature works also were of great importance for the Russian culture where 
she developed progressive and even revolutionary ideas. 
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The most important results were obtained by S. V. Kovalevskaya in the theory 
of differential equations and in rigid body dynamics. S. V. Kovalevskaya stated and 
proved the existence and uniqueness theorem for solutions of partial differential 
equations [3] (Cauchy-Kovalevskaya theorem). She also found a new method of 
integrabihty analysis (the so-called Painleve-Kovalevskaya test). The substance of 
this test is as follows. 

Let a system of differential equations 

Zi = Mzi, . . . , z.^), l<i<n (1) 

be invariant with respect to the similarity transformations 

t — > t/a, Zi — > a^'Zj, 1 < i < n. 

Let us construct a matrix (the Kovalevskaya matrix) 

Qf. 

where the constants q satisfy the algebraic system of equations /i(ci, . . . , c„) = 
—giCi (1 < ^ < JT-)- Kovalevskaya required the conditions of integrality and non- 
negativity of eigenvalues of K for the uniqueness of a general solution of (1) on a 
complex time plane. Lyapunov [27] continued these investigations of Kovalevskaya 
and related the meromorphy of a general solution with properties of the variational 
equations system. The Lyapunov method was further developed by S. L. Ziglin, who 
proved strict results concerning nonexistence of additional single- valued integrals. 

The works on rotation of a heavy rigid body about a fixed point are the most 
famous [2] (1888). In these works a new integrable case of the Euler-Poisson equa- 
tions was found and its solution in quadratures was obtained. The French Academy 
of Sciences highly appreciated this S.V. Kovalevskaya result and awarded her with 
the Prix Bordin in 1888. Let us note that the Academy of Sciences had twice an- 
nounced the prize contest on investigations on this subject, however the prize was 
not awarded. S. V. Kovalevskaya was awarded with a prize of the Sweden Academy 
of Sciences for the second memoir on the problem of rotation of a rigid body in 
spring 1889. A lot of questions concerning the motion of the Kovalevskaya top are 
still open. 



All integrable cases under consideration are well known in rigid body dynamics. 
However, not for all of them the generahzations from the algebra e(3) to the bundle 
of brackets Cx 

{Mi,Mj} = -EijkMk, {Mi.^j} = -Eijklk, {ihlj} = -xSijkMk, (2) 
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have been investigated yet. Let us note that possesses two quadratic Casimir 
functions 

c = x(M,M) + (7,7). ^ ' 

These problems naturally appear for a system which describes dynamics of a 
rigid body in an ideal fluid, connected bodies, motion of a body with eUiptic hollow 
filled in with a liquid, etc. Generalizations of known integrable cases to the algebra 
so(4) and separating variables for several of them can be found in [12] . 



1 The Kovalevskaya Top and its Generalizations 

1.1 Integrals of the Kovalevskaya Case and its Generaliza- 
tions 

S. V. Kovalevskaya found a new general integrable case of the Euler-Poisson equa- 
tions in 1888 [2] basing on pure mathematical considerations. She was developing 
ideas of K. Weierstrass, P. Painleve and H. Poincare concern the investigation of an- 
alytical extension of solutions of a system of ordinary differential equations to the 
complex time plane. Kovalevskaya assumed that a general solution does not pos- 
sess other singularities except for poles in the complex plane in integrable cases. 
That enabled to find existence conditions for the additional integral. Moreover, 
S. V. Kovalevskaya found a quite nonobvious system of variables. The motion equa- 
tions expressed in these variables are of the Abel-Jacobi type. She also obtained an 
explicit solution in theta-functions. 

In the Kovalevskaya case the body is symmetric: I\ — I2 — 2/3, where /j {% — 
1, 2, 3) are principal central moments of inertia, the centre of mass is situated in the 
equatorial plane of the inertia ellipsoid. The Hamiltonian has the form 

i^=^(M2 + M| + 2M32)+7i. (4) 

The additional integral found by Kovalevskaya can be represented as follows: 

- (^^l^ - 7i) ' + (MM, - (5) 

(without loss of generality it is assumed that the radius-vector of the centre of mass 
is r = (1, 0, 0), weight equals unit). Let us note that the integral F3 does not have 
an evident symmetry origin. 

The Kovalevskaya case (on e(3)) can be generahzed to the bundle of brackets 

(1). 
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Thus, the additional integral (the Kovalevskaya integral) will have the form 



1.2 _ l2 , L.2 

rv — A/i rvn ^ 

(6) 

ki = Mf - M| - 271 + x, 2M1M2 - 272. 

Generalizations of integrable cases can be associated with introduction of pre- 
serving integrability additional terms in the Hamiltonian without the change of 
structure of the Poisson brackets algebra (for example, introduction of a gyrostat 
or singular terms). Generalizations can be also associated with the change of the 
algebra structure without change of the Hamiltonian (for example, generalization 
to the bundle of brackets Cx), with a simultaneous change of the Hamiltonian and 
Poisson brackets, with a generalization to the case of several force fields. 

The motion equations can be generalized if one introduces a constant gyrostatic 
moment induced by a balanced rotor which is fixed in a rigid body and rotates 
with a constant angular velocity. An analogous moment occurs in motion of a rigid 
body with multiply connected (that admits a possibility of appearance of nonzero 
circulation) hollows containing an ideal incompressible liquid [18]. In this case the 
algebra and motion equations do not change. However, a linear in momenta term 
appears in the Hamiltonian. 

A gyrostat in the Kovalevskaya top was introduced by Yehia (1987) and Komarov 
(1987): 



1 (^M^ + M| + 2(m3 - + ri7i, 

F = (M^ - - 2ri7i)2 + (2M1M2 - 2ri72)'+ 

+ 4A(M3 - X){M^ + Ml) - 8riAMi73. 

In [22] D. N. Goryachev suggested a generahzation of the Kovalevskaya case at 
the zero area constant (M, 7) — 

H = 1{M^ + Ml + 2Ml) + 4 + 26i7i72 + ^2(72 - 7i) + n7i + ^^272, 



2v . . 

F = 4 ( M1M2 - 2a^ + 6173" - ^172 - r27i ) + 



7i72 , , 2 



+ Ml - 2a^^i^^ - 26273' - 2ri7i + 2r272 j , 

a, 61, 62, Ti, r2 = const. 

A physical sense of the singular term can be interpreted in quantum mechanics 
and with the help of reduction with symmetry suggested in [23] for a variant of rigid 
body motion in superposition of two homogeneous fields. Let us note that at a = 
the integrable generalization under consideration was indicated by S. A. Chaplygin 
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in his work [24], dedicated to a new integrable case of the Kirchhoff equations. 
D. N. Goryachev did not refer to [24], ahhough his work [22] appeared later. 

In the Yehia work [25] a constant gyrostatic moment along the axis of dynamical 
symmetry is added in the Hamiltonian (7) at the zero value of the area integral 

i(M2 + M| + 2(M3 + A)2) + 4 

^ 73 



H = -{Ml + Ml + 2(M3 + Xf) + — + 26i7i72 + 62(72 - 7i) + '^i7i + ^272, 
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F = 4 yMiM2 - + &173 - ^72 - ^271 j + 

+ {^Ml - Ml - 2a ^^^~^^^ - 26273' - 2ri7i + 2r272) - 
- 8A(M3 + 2\){Ml + Ml) - 16aA(M3 + 2A) {l + ^ + 

+ I6A73 (Mi(ri + 6172 - 6271) + ^2(^2 + 6i7i + ^'272)) ■ 

a, 61, 62, T2, A = const. 

In [26] Yehia suggested another generalization of the Kovalevskaya case with the 
help of the introduction of a singular (but another) term, which has not obtained a 
clear physical interpretation yet 

1 <^ 
H = -{Ml + Ml + 2M|) + ri7i + r272 + 



2 VtF+tI 

,e(M2 + M| 
W+¥ ' >i' + 7i^ 

£, ri, r2 — const. 



F = {Ml - Ml - 2ri7i + 2r272)' + 4(MiM2 - ri72 - r27i)' + +4 ^ A V + 4^—^, 



This integral curiously gives a general integrable case. 

Another generahzation of the Kovalevskaya case is associated with introduction 
of two different singular potentials in the Hamiltonian [25]: 

H = i(M2 + M| + 2M|) +ri7i + r272 + + ^1 

a, £, ri, r2 = const. 
A particular integral at the zero area constant has the form: 

F= (^M,2 - M| - 2ri7i + 2r272 - ^^^^^^P^y + 
+ 4 (^MiMs - ri72 - r27i - ^^^^ + 



VtFT^ 7i' + 7l 7I 
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Let us note that the integral F is invahd in [25] , though it is probably connected 
with a strange misprint. 

There exist generalizations of the Kovalevskaya and Chaplygin integrable cases 
which include gyrostatic terms on the entire bundle Cx- The analogue of the area 
constant is also supposed to be zero in this case 

(M, 7) = 0. 

It is more convenient to represent the Hamiltonian as follows: 

//=-((! + xa,)Ml + (1 + xa2)Ml + 2 (M3 - -) ) 



+ 



+ a27i + «i72 + + «2)73 + ^i7i + ^272- 



(8) 



This form of the Hamiltonian at A = differs from the representation indicated 
in [12] on the Casimir function \{a\ + a2)(xiW^ + 7^). The result regarding integra- 
bility of the system for A ^ is new. 
In this case the integral has the form 

K^kl + aia2kl - X{ki, k2} - AX\M^ + M|), 

ki = aiM^ - Q.2M2 - (oi - 02)73 " 2(7iri + 72r2) + x- 



aia2 



u HI T^i ^"1^172 + "2^^271 ,0 nr2 
k2 — M1M2 — 2 h 2a;- 



aia2 C11C12 

2 2 

^^2} = AuJa^Ml - a2Ml - ^^i!l±^^ _ 4^3(^0^ _ a2){M^ii - 
\ aia2 / 

-2(Miri + M2r2)). 

(9) 

Remctrk. The motion equations for K can be represented in the form 

k = {K, H} = (ai - a2)(M, j)F{M, 7). 

Thus, this case becomes a general integrable case at ai = 02 that corresponds to the 
gyrostatic generalization of the Kovalevskaya case. 



1.2 Separation of Variables on the Bundle Cx 

Our reduction to the Abel equations on the bundle of brackets uses arguments of 
G. K. Suslov, who suggested a new method of integration of the Kovalevskaya case 
in his famous treatise [4] . 
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Let us set new variables 



^1 = Ml + 1M2, Z2 = Mi- iM2, 

Ci = ki + ik2, (2 = h- 'ik2, 

C1C2 = k'. 

Let us take advantage of the motion equations for zi , Z2 

izi = M3Z1 - 73, -iz2 = M3Z2 - 73 
and formulate 71, 72, 73, M3 as 

72 = - ^ {4 - zl) + 1 (Ci - C2) , (10) 

.Z1Z2 + Z2Z1 ,Zi+Z2 

73 = i , Ms = I . 

Zl — Z2 Zl — Z2 

Then insert the obtained expressions into the integrals (2) and the Hamiltonian 
(4). Solving the obtained equations, one finds 

il = -Rl - ^ {Zi - Z2f , 

ZI^R2-^^{ZI-Z2)\ (11) 

Z1Z2 = - R - ^{2h - x) {zi - Z2f , 

where 

R^R{zi,Z2)^\zfz^-!l{zf + zi) + l{zi + Z2) + ^-c + xh-^, 

[12) 

Rl = R{Zi,Zi), R2 = R{Z2,Z2), 

here /i is a constant of the energy integral (4). 

We still need to exclude d, ^2 from the obtained equations with the help of the 
Kovalevskaya integral (6) 

{Rl - zl) {R2 - zl) = ^-^{Zi - Z2)' = ^^{Zi - Z2)\ 
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Rearrange terms in the obtained relation and represent the latter in the form 

{ ■^l _|_ _£2_^ ^ — { ZlZ2 _|_ 1 V _ fc^(zi-22)'* _ f 

2 2 ^^^^ 

( _zi Z2_\ _ ( ziZ2 _ 1 ^ _ k^izi-Z2)^ _ [■ 

VVRT y-^/RiB^ ) I6R1R2 •'2' 

where Z1Z2 are substituted from (11). 

The last step is the introduction of the Kovalevskaya variables by formulae 



R — \/ R1R2 R + \/-Rii?2 



Let us express R = {si + S2){zi — z^f ^ \/R\R2 = (s2 — — -22)^ from these 

relations and insert them in the right-hand sides of (13). We obtain 

f _ f{si) f _ f(s2) 

f{s) = {2s + l{2h-x) + \k) {2s + 1(2/1 - x) - \k) . 

The relations for the left-hand sides of (13) define a transition to curvilinear 
coordinates 

dzi dz2 dsi dzi dz2 ds2 



where (p{s) is a polynomial of the third degree, 

^{s) = 4.=^ + 2hs' +(^^i2h-xr-^ + ^y + ^. 
Inserting (16) into (13), we obtain the motion equations in the variables si, S2 

Si — , S2 — . [i-l) 

Si — S2 S2 — Si 

The polynomial ip{s) can be obtained with the help of standard methods using reduction of 
elliptic integrals to the standard from. Let us describe the Suslov method [4] here. 
The Kovalevskaya variables are solutions of the quadratic equation 

Q{zi, Z2, s) = {zi - Z2fs^ -Rs + G = 0, (18) 

where 

G = = -^^1^2(^1 + ^2) + {{2h - xf -k^+ 4c) {zi + z^f - ^(zi + Z2) + J. 
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Let us calculate quadrates of derivatives , (^) ' ^^'^ using (18) exclude s,zi,Z2 

from them respectively. By a direct calculation we verify that 



Let us construct an exact differential of Q 

aQ = -T^dzi + ——dz2 + -r^ds = 0. 

Dividing by ^J Lp{s)RiR2 and taking into account a possibility of rooting with different signs, we 
obtain the relation (16). 

1.3 The Haine-Horozov Transformation for the Kovalevskaya 
System 

Let us consider another method of integration of the Kovalevskaya system. This 
method uses a coordinate transformation analogous to the one indicated in [5] . The 
Kovalevskaya system reduces to the Neumann problem on a motion of a point on 
S"^ with the help of this transformation. That allows to write down a Lax pair for 
the Kovalevskaya top and implies a relation between integrals of both systems. 
Let us denote the motion integrals of the Kovalevskaya system on Cx as follows: 

2H = ]-{Ml + M| + 2M|) + 271 = h, 
ci^{M,^), C2^x{M, M) + {j,j)-e, (19) 

k = 71 + x] + — 72 



4 

and transform variables by formulae 

_ , _ 4Mi73 -M3 {Mf +M| -4) 

~ 4M2 ' ^ ~ 8M2 ' 

- -^t' ^2 = (20) 

_ , Mf+M|-4 , _ , 4Mi73-M3(M-P+M|+4) 

Then after the coordinate transformation (20) the motion equations can be writ- 
ten as 

xix — h] 

p = lxp, l = Qpxp + {pi + ip3)y, (21) 

where the matrix Q has the form 

^ I C2-I -ici i{l + C2 
Q = - I -ici -h ci I , (22) 

i(l + C2) Ci 1 - C2 



4 

h, Ci, C2, are determined by (19), the vector y = {—ip2, i{pi + ips), P2) 
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The motion equations (21) can be regarded as the system on e(3) with the 
Hamiltonian 



{I, I) + (Qp, p) + ^ ^ \ pi + ipsf - ^ = -4 > 

i. e. as the Neumann system. Thus, the Kovalevskaya system on the bundle (1) can 
be integrated by reduction to the Neumann system. The integration of the latter is 
performed by standard methods. 

Let us note that the Kovalevskaya system can not be reduced to the Neumann 
system with the use of the transformation (20) after introduction of a gyrostat. 



1.4 A Generalization of the Kolosoff Analogy 

G.V. Kolosoff considered a transformation in [7] which involves coordinates and 
time and reduces the Kovalevskaya problem on e(3) to dynamics of a point on the 
Euchdean plane in in a potential field, for which variables can be separated. 
This transformation is used to introduce action-angle variables in the Kovalevskaya 
problem in [8]. 

Let us consider an analogous procedure for the Kovalevskaya problem on the 
bundle (1). In this case an analog of the Kolosoff transformation gives dynamics of 
a particle on an axisymmetric nonconstant curvature surface. 

Following [8] , let us express the Hamiltonian from the motion equations. Let us 
change coordinates in (17) 

s,^Si-^{h+^), 1=1,2 (23) 
and represent them in the form 

M) ' 

/(«)=4(.-| + ifc)(5-|-lfc), 

g{s) = As^ - {2h + lx)s^ + {c-k'' + + x, (24) 

X = -{k^h + 2/2 - ch) + ^{h^ + e-c)- (^y. 

Remark. The form of the coordinate change (23) is dictated by the requirement 
that the energy constant in g{s) was contained in even powers of s only. 
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Substraiting the first equation from tlie second in (24), we eliminate the constant 
X and express the energy from the obtained relation 



, 2(5? + 5i., + 5i) + |(c-P + ^) ^ ^ 



Si + S2 

After the change of time dr = 2(si + S2) dt and transition to the canonical 
momenta pj = f7)'5i = ^) ^ = 1)2, one obtains a system with separable variables. 

Let us consider the variables si — |, S2 — f as elliptic coordinates on a plane 
with the following relation with the Cartesian coordinates v 

2 / x\ f x\ k 

We obtain a system on a plane with the energy determined by the relations 

H =T + U, 

^ _4p'^ - 2uk + c + 3xp + _ 2(p2 + pI)-P + c + 3xp + x^ ^^^^ 



( x\ ( x\ k 
(^^-4) (^^-4) + 2' 

2 IT. Fx /A;2 



2p + x 2p + 2; 

p —\f\TTv^ 1 p\ — \l {u — ky + v^. 

The system (26) describes the motion of a point along a curvilinear surface in 
a potential field. The Gauss curvature can be evaluated by the metrics which is 
defined by the kinetic energy T 

K- 



(2p + a;)3" 

1.5 The Act ion- Angle Variables for the Kovalevskaya Top 
on the Bundle 

The method of introduction of the action-angle variables developed in this paper 
is analogous to the method considered in [8]. In comparison with a usually cited 
procedure due to A. P. Veselov and S. P. Novikov, this algorithm is more natural and 
uses an ordinary method of introduction of action-angle variables for the systems 
with separating variables. 
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Our algorithm consists of the following steps: 

1. Construction of the Abel variables Si, S2, that commute in the original Poisson 
structure: {si, S2} = 0. 

Remark 1. The existence of commuting set of the Abel variables si, S2 can be 
established with the help of arguments due to [17]. These arguments are connected 
with the reduction of the equations (17) to the standard form on a torus. 

2. Introduction of the canonical momenta pi with the help of the energy equation 
in the variables Sj, Sj. These momenta must satisfy an additional requirement that 
the system pj, Si should be separable. 

3. Having a set of separated variables, the action-angle variables can be intro- 
duced in accordance with a well known algorithm. 

It is easy to verify that the Kovalevskaya variables (14) si, S2 commute. They 
also satisfy the equations (17) 

. _ \/ f{si)ip{si) . _ v//(s2)'/?(g2) 
si — , S2 — , y^t) 

Sl - S2 Si- S2 

where 

fis) = (^2s + ^(2/1 -x) + J/c) (^2s + ^{2h -x)- \k ^ 

.3 2 n,o. .2 c\ P ' ^^^^ 

cp{s) = As' + 2hs' + i^-{2h -xr-- + -js + -. 

(The variable Si varies from to 00, and S2 parameterizes a circle.) 
Let us extract a motion integral >i = jq {{2h — x)"^ — k"^) in (28): 

f{s) = 4s^ + {2h - x)s + X, 

c f (29) 
93(5) = s(4s^ -h 2/is + X + -) + — 

and exclude the variable x from (29). The obtained equation express the energy h 
as a function of s, and s, : 



/I - -2(^Si + S2j + + 4 + , .„„x 

_ 16sfa;+4s,c+;^ _ (st-S2)si ^'^^ > 

64si ' - 2^i ■ 

Let us introduce the conjugated momenta Pi instead of the velocities sf 

where F[si) is an arbitrary function of Sj. Note that the addition of F[si) does not 
change the motion equations (the canonical transformation). 
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We obtain 

p,^^lnfi±^^. (32) 

Let us write down the Hamiltonian of the Kovalevskaya system (30) in the vari- 
ables Si, Pi 

/2 aiC0sh(2pi^) - a2COs(2p2y^) 
-si- S2 + h-H . (33) 

8SiS2 8 Si - S2 

The variables Si are separating for the Hamiltonian (33). Introducing the con- 
stant of separation >ci, we obtain two equations, which can be integrated separately 
from each other: 

2sj + si (/i - f ) + + xi = ai cosh(2piv^), 
2s| + S2 (/i - f ) + + xi = 02 cos(2p2^=s^). 

Inserting (32) in (34), we obtain the Kovalevskaya equations (27). The separation 
constant Xi is connected with the constant x by the formula x = 2xi — |. 
Thus, the action in the variables s equals 



1 f 



2n 



f 



p{s) ds, 



i. e. is proportional to the area, restricted by a phase curve. The integration region 
depends on the parameters and values of integrals of the system. 



2 The Goryachev-Chaplygin Case (1903) 
2.1 Integrals and Separation of Variables 

The body is dynamically symmetric in this case: I1 — I2 — 4/3, the center of mass 
is situated in the equatorial plane of the inertia ellipsoid. The Hamiltonian and the 
additional integral have the form 

H=]^{Ml + Ml + AMl)+^,, ^^^^ 

The gyrostat in the Goryachev-Chaplygin case was introduced by L. N. Sretensky 
(1963): 

H=^ (^M^ + Ml + 4(m3 - + ri7i, 

F = (2M3 - k){M^ + Ml) - 2riMi73. 

13 



Let us note that the gyrostatic moment is directed along the axis of dynamical 
symmetry in generalizations of the Kovalevskaya and Goryachev-Chaplygin cases. 
A separation of variables in the Sretensky case (a generalization of the Goryachev- 
Chaplygin case) is described in [19]. 

A generalization of the Goryachev-Chaplygin case at the zero area constant was 
suggested in the work of D. N. Goryachev [20]. The generahzed Hamiltonian and the 
integral have the form 



H = \{Ml + Ml + AMI) - ni + 4' 



F = Ms ( + Ml + 2^ ) + nsMi, 



a 

a, r — const. 



(36) 



I. V. Komarov and V. B. Kuznetsov in [21] added a constant gyrostatic moment 
to this Hamiltonian (also at {M, 7) = 0) analogous to the Sretensky generahzation 



H = 1{M^ + Ml + 4M|) - r7i + AM3 + ^, 



F = ( M3 + ^) + ^2 + + ^73Mi, (37) 

a, r, X = const. 



We shall give an analog of the Andoyer-Deprit variables for the bundle of brackets 
Cx. Let us use of the following sequent of embeddings M} C 5*0(3) C Cx- 
Let Ms be equal to one of the momentum variables 

M3 = L. (38) 

Its canonically conjugated coordinate /({/, L} = 1) on so(3) can be found by inte- 
gration of the Hamiltonian flow with the Hamiltonian Ti — L 

^^{M„L}^M,, ^^{M,,L}^-M„ ^ = {M3,L} = 0. (39) 

Further taking into account the commutation relation {M2, M2} = — M3 we obtain 



Ml = VG^-L^sinl, M2 = VG^-L^ cos/, (40) 

where — Mf + M| -|- M| is a Casimir function of the subalgebra so(3). Let 
us choose G as the other momentum variable. Then the corresponding flow on the 
total algebra has the form 

dM ^ dp 1 ,,,, 
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where g is the variable canonically conjugated to G. 

Remark. The choice of G as a new canonical variable (not G"^) is determined 
by the fact that the corresponding variable g varies in the range from to 27r and 
does not depend on the value ofG. 

Accordingly to (41), M does not depend on g. Using (41) for p and the Casimir 
functions 

(M, p) ^H, p'^ + xM^ = c 

we obtain 

H a 

P= —^M+-{M X easing + G'M x (M x eg) cost/), 

c - xG^ - ^ ^^^^ 
= — ^3 = (0,0,1). 

Thus, (38), (40), (42) define symplectic coordinates on the entire bundle C^- 
These coordinates at x = pass into known Andoyer-Deprit variables in rigid body 
dynamics. 

We shall find a generalization of the Goryachev-Chaplygin case on using (38), (40), (42). 
Let us take the Hamiltonian in the form 

l-{G'^ + 3L'^) + XL + a{coslcosg + ^smlsmg), (43) 

where a, X are constants. 

In comparison with [17] , a linear with respect to L term describing the gyrostat 
is added in (43). 

The system (43) admits a separation of variables. Indeed, let us perform a 
canonical change of variables 

L = pi+p2, G = pi-p2, qi = l + g, q2 = l-9- (44) 
The Hamiltonian (43) can be represented as 

Ipl-Pl M-PI. ^ ( ■ , - A (AV,\ 

ri^- A \ (pismgi +p2smg2)- (45) 

2pi-p2 P1-P2 P1-P2 

Let us express the Hamiltonian (43) via the variables M, p at the zero area constant 

(M, p)^H^O. 

We obtain 

n = -(M^ + M| + 4M|) - AM3 + aP-. 
2 \P\ 

The additional integral in this case has the form 

n^{Ms- ^){M^ + Ml) - aMi^ 
2 \P\ 

For the algebra e(3) \p\ = 1, thus we obtain the Goryachev-Chaplygin case. 
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2.2 The Action- Angle Variables for the Goryachev-Chaplygin 
Case 

The variables qi, q2 (44) are already separable. Moreover, they commute with each 
other as it is easy to verify. 

Let X denote the separation constant. Then 



2 

pi 



,3 

— Xpf + api sin qi — hpi — x, 



Pi 



2 Xpl - ap2 sin q2 - hp2 = x. 



The action in the variables s equals 



I = p{s) ds. 



(46) 



Remcirk 2. Since 

dH api cos qi 

Pi = -jr~ = ' (47 

oqi Pi - P2 

using (46) we obtain the Abel equations in the form 

2{pi-P2y (48) 
$(^) = Aa'^z'^ - (2x -z^ + 2Xz'^ + 2hzf. 



3 The Chaplygin Case 

3.1 The Integrals and the Separation of Variables 

This case discovered by Chaplygin in 1902 [24] is a particular integrable case at the 
zero area constant (iVf , 7) = and a quartic integral. This system is similar to the 
Kovalevskaya case in the Euler-Poisson equations. 

The Hamiltonian and the additional integral of the system have the form [24] 

//=^(M,^ + M| + 2M|) + ^(7,^-7D, 

F = (Ml - Ml + -7^)2 + AMlMj. 
a 

Let us consider the integration of the Kirchhoff equations using the method 
of separation of variables in the Chaplygin case on the bundle of brackets Cx 
(the integration in this case has not been indicated before). The separation of 
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variables in this problem on was performed by S. A. Chaplygin [24] e(3), and by 
O. I. Bogoyavlensky [12] on so(4). 

The Hamiltonian and the integrals on the bundle have the form 

H = KaaMf + q;iM| + (ai + q;2)M| - (oi - a2)(7? - tD) = |, 
J2 = X ( + M| + M|) + 72 + 7I + 72 , 
J3 = Mi7i + M272 + M373 = 0, 
J4 = (aiMf - a2M^ - (oi - 02)7!)^ + 4.aia2M^M^ = fc^, 

where ai = 1 — xai, ^2 = 1 — 2:02 (A; > 0). Let us introduce the separating variables 
si and S2 by formulae [12, 24] 

„ u+k „ u—k 

(50) 

u = aiMf + 0:2^^25 V = {a2 — 01)73- 
The evolution of si and S2 is determined by degenerate Abel-Jacobi equations: 

h = -^J{l-sl){6l-P,S,), S2 = -^{l-sl){S2-P2S2), (51) 

where 

Si — 'mi{h + k) — m2 J2, S2 — mi{h — k) — m2 J2, 
(3i = x(ai - a2){h -\-k) - m^J2, P2 = x(ai - 02) (/i - k) - 1713 J2, 
mi^ai + a2, m2 = x(ai - 02)^, 777.3 = (oi - 02) (cti + 0:2). 
Thus, the motion equations can be integrated by means of elliptic functions of time. 

3.2 The Action-Angle Variables for the Chaplygin System 
on the Bundle 
of Brackets 

The variables si, S2 (50) commute at the zero value of the integral (Af , 7) = 0. 
We find the energy E — 2h as a, function of si, S2, s'l, S2 from (51): 



^ £1^ ^ £2^ 2J2m3 ^ 

(1 - sf)(siai2 + mi) (1 - s|) (52012 + mi) ai2 
^ J2(m2ai2 + mim3) / 1 ^ 1 



(52) 



ai2 \siai2 + mi S2O12 + mi 

where ai2 = x{a2 — ai). Introducing the conjugated momenta (31) 

Pi 



[1 - s'^){siai2 + mi) 
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we obtain the Hamiltonian in separated variables 

if = |(1 - Si)(siai2 + mi)pl + |(1 - sl){s2ai2 + rni)pl - 'H^I^^j^ 

^ J2(m2ai2 + mims) / 1 _^ 1 ^ ^^"^^ 



ai2 \Siai2 + mi 52012 + "7.1 

As it may be expected, the variables are separated in the Hamiltonian (53): 

-(1 - sj(siai2 + mi)pi H ■ — = x, 

4 ai2(siai2 + mi) 

1 , X 2 , ^2(^2012+^17713) 2727773 

-(1 - S2) (-52012 + ^ii)P2 H 7 : ^ = ^ H ^■ 

4 012(^2012 + "7i) ai2 

The action in s equals 

I = p{s) ds. 

4 The Bogoyavlensky System 

The particular Bogoyavlensky case [12] with the Hamiltonian 

H^^ (^aiM^ + a2M^ + 03^^32 + ^((02 + 03)7? + (oi + 03)72 + (oi + 02)7!)) = ^ 

can be integrated by means of elliptic functions. The system possesses the integrals 
J2 = x(Mf + M| + Ml) + 72 + 7I + 7I, 
J3 = Mi7i + M272 + M373 = 0, 

■Ji = ((03 - 02)7? + (oi - 03)7! + (oi - 02)7!)^ + 4(03 - ai)(o3 - 02)7?7| = 

On so(4) the integration method is developed in [12]. If we introduce the variables 

u + h u — h 
si = , S2 = , 

V V 

where u = {a^ — 02)7^ + (03 — 01)7!, = (oi ~ 02)7!, motion equations can be 
reduced to the form 

s^^-^{l-sl){Si-Pis,)/2, s2^-^{l-sl){S2-P2S2)/2, (55) 

and the constants 5i, (5i are expressed as follows: 

5i = {Ji + h/2)ni + 72772, 62 = {Ji - h/2)ni + J2772, 
/5i = (oi - 02)(-^2a3 - Ji- h/2), p2 = (fli - a2){J2a3 - Ji + h/2), 
Til = 2a3 — tti — 02, 772 = 2aia2 — 03(01 + 02). 
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Remark 3. Let us note that (51, 55) are degenerate Abel-Jacobi equations, i. e. 
each of them depends on a unique variable Si or S2 only, and a two-dimensional 
Abel torus splits into one- dimensional tori. 



5 Integrable Systems on a Sphere with a Cubic 
Integral 

(the Gaffet System) 

In the B. Gaffet work [13] an integrable case on the two-dimensional sphere 7^ + 
7I + 73 = 1 with the potential U = 0(717273)"^^^ is found. 

This integrable system on e(3) and at the zero value of the area integral possesses 
the integrals 



Let us note that the integrability of the system and the form of integrals (56) 
are preserved on the total bundle of brackets and at the zero value {M, 7) = 0. 

Despite some particular results [14, 13], an explicit integration of (56) has not 
been performed yet. L — A pair for this system is indicated in [15] and has the form: 
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where y 
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